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Renormalization constants of local operators within the Schrodinger 
functional scheme 

A. Shindler a * 

a Dipartimcnto di Fisica, Universita di Roma Tor Vergata 

We define, within the Schrodinger functional (SF) scheme , the matrix elements of the twist-2 operators 
corresponding to the first two moments of non-singlet parton density, and the first moment of singlet parton 
densities. We perform a lattice one-loop calculation that fixes the relation between the SF scheme and other 
common schemes and shows the main source of lattice artefacts. Few remarks on the improvement case are 
added. 



1. Introduction 

The accurate knowledge of hadron parton den- 
sities is an essential ingredients for the experi- 
mental text of QCD at the accelerator energies. 
Their normalization is usually obtained from a fit 
to a set of reference experiments and used for pre- 
dicting the behaviour of hard hadron processes 
in different energy regimes. The calculation of 
the normalization needs non-perturbative meth- 
ods. These computation, expecially for the higher 
moments, [jj] can reduce, for example, the uncer- 
taintes on the gluon parton densities at values of 
the Bjorken x larger than 0.5. These calculation 
are made mainly by two groups (see for example 

It's well known that, to match the scheme of 
the non-perturbative simulation and the scheme 
where the experiment made the comparison with 
theory, it is necessary a lattice perturbative com- 
putation of the renormalization constants of the 
operator in the same scheme where the operator 
is numerically computed. 

The perturbative and non-perturbative || 
calculation are made in the Schrodinger func- 
tional scheme (SF). The Schrodinger functional 
has been discussed extensively in the literature 
(see || for reviews). Among the advantages of the 
method, we only quote the possibility of perform- 
ing the computations at zero physical quark mass 
and of using non-local gauge invariant sources 
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for fermions and gluons without need of a gauge- 
fixing procedure. In our particular case, we ex- 
ploit both features. 

2. Non-singlet 

We define the SF correlation function of the 
first and second moment of the non-singlet parton 
densities by the observables: 

/o 12 (x ; p) = / 2 (x ; p) = -a 6 ^ * 

(^Wth B 2] irV(^)C(y)rir 3 c(z)) 

/o 123 fro; P) = fs(x ; p) = -a 6 ]T x 

(^(x) 7[1 S 2 S 3 ] irV(x)C(y)irr 3 C(z)) 

where the contraction of the classical fields is non- 
vanishing if the matrix T satisfies: rP T = P±T, 
where P± = h(l ±70) and p is the momentum of 
the classical field sitting on the boundary. 

The quantities £ are the response to a variation 
of the classical Fermi field configurations on the 
boundaries. We take the limit of massless quarks, 
but some care should be taken so as to ensure this 
limit at order g 2 . The breaking of chiral simmetry 
of the Wilson action entails a non-zero shift of the 
quark mass from the naive value at order g 1 . 

The matrix element of the operator for the first 
moment involves two directions and three for the 
second moment. These directions must be pro- 
vided by external vectors: we have chosen to ob- 
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tain one of them from the contraction matrix T, 
i.e. from the polarization of the vector classical 
state r = 72, and the remaining ones from the 
momentum p of the classical Fermi field at the 
boundary. 

To compute the renormalization constants of 
the operators is necessary to remove the renor- 
malization constant of the classical boundary 
sources £. Following ref. ||, this is represented 
by the quantity called fx . Both fa and fx are nor- 
malized by their tree-level expressions. We define 
the renormalization constants such that operator 
matrix elements (here briefly indicated with O) is 
equal to its tree- level value at fx = 1/L: 

O r {h) = Z{an)- l O haxe {a/L) (2) 

with Z(a/L) defined by 

O barc (a/L) = Z(a/L)O trcc (3) 

At the one loop the observable can be 
parametrized as: 

Z(pL, x /L, a/L) = 1 + g 2 Z^ (a/L) (4) 
with 

Z {1 \a/L) = 

. , , /n , k b k + c k \n(a/L) (5) 
b a + c \n{a/L) + ^ a K \ 

k=X 

The study of the scale dependence of the renor- 
malization constant is then equivalent to that of 
the dependence upon the lattice size L, provided 
the external variables, upon which the matrix el- 
ements depend, scale like the basic lenght L. We 
made the calculation with the minimum value of 
the momentum pL = 2tt, and with a "finite size" 
momentum 8 = 0.1/ L, which escape the momen- 
tum quantization rule. Both the computation are 
made with Xq = L/4, and xo = L/2. . The fit 
(the fitting procedure is described in ref. [Q) to 
the N = L/a dependence of our results has been 
made by the expression 

Z^(N) = 

* + fttaW+ £*±£5HQ (6) 

k=X,2 



Table 1 

Values of the constants of the operators renor- 
malized in the two definitions with the real mo- 
mentum p or the "finite size momentum" 9 dif- 

ferent from zero. 

Moment Definition Constant 
First x = L/4(p ^ 0) B = 0.2635(10) 
First xo = L/2(p ^0) B = 0.2762(5) 
Second xo = L/i(p ^0) B = 0.1875(20) 
Second x = £/2(p^0) Bp = 0.1895(50) 
First xo = L/4(9 + 0) B = 0.12180(15) 
First xo = L/2(6 ^0) B = -0.23340(15) 
Second xo = L/4(6 ^ 0) B = -0.06080(55) 
Second xp = L/2(9 ^0) B = -0.5675(20) 



excluding higher-order terms. Tabic |l] contains a 
summary of the constants Bo for the operators, 
after removing the external legs renormalization, 
in the various cases that we have discussed. The 
lattice result converge much faster in their contin- 
uum limit, using a "finite size" momentum, con- 
firming the momentum quantization as a major 
source of lattice artifacts (see ref.0]). We made 
the same computation with clover action J?J . Also 
in this case lattice artefacts start at order a for 
both coefficients, because the 0(a)-improvement 
of the action should be complemented with the 
improvement of the operators and of the bound- 
ary counterterms in order to lead to a full cancel- 
lation of effects appearing linear in a. The Feyn- 
man rules for this computation can be easily de- 
rived using ref. || . The calculation is done to the 
order c 2 sw . It's important to stress out that to ex- 
tract the finite constant to this order we must do 
the computation to the same order of fx and of 
the mass shift. The fitting procedure is the same 
that in the Wilson case. The finite part of the 
renormalization constant now will be defined as 

B =B^+B^c sw +B^c 2 sw (7) 

Doing the subtractions of fx and of the mass shift 
we obtain the following results: 

B [ o ] = -0.0327(6) B^ ] = -0.0327012 

B {2) = -0.005725(9) B {2) = -0.005726859 
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Where Bq , and Bq are the coefficients com- f gg (x ) — a 6 y^(Tr{Fi p F P 2}Gi2(y, z)) 



puted in a standard lattice It's clear that the 
parts of the finite constant proportional to c sw 
and to c 2 sw are the same, within the errors, in the 
two schemes. So the difference between the total 
finite constants in the two different schemes will 
be indipcndcnt on the fermion action used. 

3. Singlet 

We can reproduce the same calculation in the 
flavour singlet sector. In the singlet case the mix- 
ing under renormalization of the operator force to 
compute 4 correlation functions, involving gluonic 
and fcrmionic boundary state. The 2 SF correla- 
tion function involving fermionic boundary state 
are: f qq {xo; p) which is essentialy the same of / 2 , 
and 
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fgq(xo;p) = 

-a 6 J2 c ip(y " z) (Tr{F 1 ,F, 2 }C(y)rC(z)) 



(9) 



To study the other two correlation functions it is 
necessary to fix the gluonic boundary field. We 
use again the particular gauge group that leaves 
invariant the Schrodinger functional. We define 
the Big Teeth state as follows 



T ■ 



A(y) = { n u °to + fca o)}^(y + t°) 



k=0 



(10) 



{ II U^iy + ai + kaO)} 

MM 

We define then the gluonic boundary state as 
G 12 (y,z) =Tr[ 5l (y) 52 (z)] (11) 
where 

1 



i(y) 



2iago 



(A(y)-A(y)- 1 } 



(12) 



The 2 remaining correlation function are defined 
as: 

fqg(xo) = 

- a e £<^(z)7[i B 2] i>(x)G 12 (y, z)) ^ 13 ) 



It's important to note that the SF scheme allow 
us to construct a gauge invariant gluonic source 
projected at momentum with a tree-level well 
defined in the continuum limit. This is achieved 
by the two different polarization of the gluonic 
source. This improve the convergence to the con- 
tinuum limit as it's been remarked in the non- 
singlet calculation. The results for f gq and f qg 
are the following: 

B$ = 0.0261(3) B q °J = 0.02919(1) (15) 

It's currently under way the computation of f gg 
and of the gluonic source renormalization. 
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